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ABSTRACT 

In this Letter we model the chemistry of DCO + in protoplanetary disks. We find that the over¬ 
all distribution of the DCO + abundance is qualitatively similar to that of CO but is dominated 
by thin layer located at the inner disk surface. To understand its distribution, we investigate 
the different key gas-phase deuteration pathways that can lead to the formation of DCO + . Our 
analysis shows that the recent update in the exothermicity of the reaction involving CH 2 D + as 
a parent molecule of DCO + favors deuterium fractionation in warmer conditions. As a result 
the formation of DCO + is enhanced in the inner warm surface layers of the disk where X-ray 
ionization occurs. Our analysis points out that DCO + is not a reliable tracer of the CO snow 
line as previously suggested. We thus predict that DCO + is a tracer of active deuterium and in 
particular X-ray ionization of the inner disk. 

Subject headings: protoplanetary disks — astrochemistry — ISM: abundances — stars: formation 


1. Introduction 

The D/H ratio is a key tracer of history and 
origin of chemical and physical processes that 
occurred in the early phases of protoplanetary 
disks, in which plan ets m ay form. In that light, 
van Dishocck ct al. (2003) have suggested that 
deuterated chemistry pathways are active during 
the early phase of planet formation which may al¬ 
ter the deuterium enrichment supplied by the ear¬ 
lier, pre-stellar stages. These deuterium enhance¬ 


ments may be provided to ice-rich planetesimals 
during formation. Thus understanding deuterium 
fractionation is key to unravelling the clues of the 
D/H ratios in t he various bodies found in our solar 
system (e.g. iMumma fc Charnlcvll20li l. 

In the physical conditions (temperature, den¬ 
sity, ionization) prevailing in circumstellar disks 
surrounding T Tauri stars, deuterium fractiona¬ 
tion is expected t o be dominated by ion -molecule 
reactions (see e.g. Aikawa fc Herbst 1999l l. Among 
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deuterated ions, DCO + has been detected in var¬ 
ious disk environment such as toward the pro 
topla n etary disk o f DM Tau dGuilloteau et al 


and 


2006; Oberg et al. 2010l 2011a : Teague et al 


20151) . th e disk around the classical T Tauri star 
TW Hya (Ivan Dishoeck et al.ll2003 ; Qi et al.ll200S ; 


Oberg et al. 2012) as well as in the disk surround¬ 


ing HD 163296 (jMathcws et al. 12013 ) but also 
in the following other disks, LkCa 15, IM Lup, 
AS 09, and V404 6 Sgr ( SMA DISCS project, see 
Oberg et alJl2010L l2011bh . 

The key gas-phase deuterium pathways that 
are active at cold tempera ture (typically T<50K) 
are th e following reactions (iVV'atsorJI 1 97fit IWootten 
19871) : 


h t- 


HD ^ H 2 D+ + H 2 + 232K 


and then, 


H 2 D+ 


CO -+ DCO h 


h 2 


(1) 


( 2 ) 


Reaction [T] results in H 2 D + enhancement as it 
is driven mainl y forward ( i. e. lef t-to-right) at low 
temperature ( Pagani et al.lfl992 fl. Then by react¬ 
ing with CO, more specifically by transferring a D- 
atom to CO, the H 2 D + ion can lead to the forma¬ 
tion of DCO + (reaction. [2]). Due to the low tem- 
peratu re condition that is fav orable to reactions [l] 
and l2l [Mathews et al.1 ( 2013 ) have suggested that 
DCO + is a good tracer of the disk midplane near 
the CO snow-line. The latter represents a line 
dividing regions where the disk is warm enough 
for CO to be in the gas phase and regions where 
the temperatures are low enough for CO to be de¬ 
pleted on the grains. 

However, modeling the chemical structure of 
a protoplanetary disk is not trivial as the chem¬ 
istry depends of the physical env i ronm ent (e.g. 


iDutrev et al.l f2014i lAikawa et al.l 2002h ; and 


DCO + may form independently of H 2 D + though 
reactions with CH 2 D + unde r warmer cond itions 
(typically T> 50K, see e.g. IWootten 19871) . In 
particular, in warm conditions the CH 2 D + ion is 
mainly formed by 

CH+ + HD ^ CH 2 D+ + H 2 + AE. (3) 

Reactions involving CH 2 D + or one of its products, 
such as reaction 


CH 4 D+ + CO -+ DCO+ + CH 4 , 


(5) 


ch 2 d+ + o -> dco h 


h 2 


(4) 


can significantly contribute to the formation of 
DCO+. 

Until now, the zero point vibrational energy 
difference between the products and reactions, 
AE, of reaction p was esti mated to be 370 K 
( Smith ct al. 1982f b However, Roueff et al.l ( 2013 1 
have recently revised the energy budget of reac- 
tion[3]and estimate an exothcrmicity AE of 654 K, 
which is about a factor 2 larger than the previous 
estimate. This finding implies that CH 2 D + can 
become enhanced even at temperatures as high as 
300 K. As a consequence, the DCO + formation 
rate could also be increased. Incidentally, using 
Plateau de Bur e Interferometer obse rvations of the 
DM Tau disk, Teague et al. ( 2015ll have recently 
pointed out that DCO + can be used as a potential 
tool to diagnostic ionization in disks. This leads 
one to ask how the change of exothermicity af¬ 
fects the deuterium chemistry and hence, whether 
DCO + is a reliable tracer of ionization in disk or, 
of the CO snow line as previously suggested by 
Mathews et al. ( 2013l f. 

In this Letter, we discuss the distribution of 
the DCO + ion i n disk according the new assump¬ 
tions set bv iRoucff ct al. (2013 ). In Section [2] we 
describe the physical and chemical model used in 
this study. Results are given in Section[3]and their 
implications are discussed in Section [I] 

2. Modeling 

We use the (1+1)-dimensional disk model of 
Fogel et all ( 201lh to investigate the main forma¬ 
tion pathways of DCO + . The physical structure 
and reaction network are described below. 

2.1. Disk Structure 

In this study, a two-dimensional and az- 
imuthally symmetric disk physical structure that 
represents a T-Tauri disk has been adopted (i.e. 
structure of the TW Hya T-Tauri system see 
Cleeves et al.l 12013 : ISchwarz & Berginl 20l4 and 
references therein). Being the same settled disk 
struct ure as the one adopted by Schwarz fe Berginl 
(|2014l) . we simply refer to this paper for a complete 
description of the density/temperature structure 
and grain distribution. 
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2.2. Chemical reaction network 


3. Results 


We use the same deuterat ed chemical net¬ 
work as in lAikawa et al. ( 2012 ), who worked on 
star forming cores, except for the following up¬ 
dates and modifications needed for disk chem¬ 
istry: i) the energy barriers of the endothermic 
reactions involving th e CHo D + ion are updated 
following iRoueff et al.1 ( 20131) . ii) The X-ray ion¬ 
ization of H 2 and He along with UV photo ly¬ 


sis induced by X-rays are included (|Fogel et al 


20111) . in) Lya photodcsorption radiation on dust 


grain is taken into acco unt Isee lFogel et aljl201li 
Bethell fe Berginl 2011 . for further details), iv) 


CO, H 2 , HD and D 2 self-shie lding are include d 
in a time-dependent way (see Fogel et ahl 1201 lh . 
v) Finally, grain surface reactions are limited to 
the formation of H 2 , OH and H 2 O since gas-phase 
chemistry prevails in the DCO + formation. 

The network contains in total 1565 species and 
43353 reaction^. Among the reactions listed 
above, the model also includes grain surface re¬ 
combination reactions, adsorption, photodesorp¬ 
tion, two body reactions (including ion, molecule, 
negative ion, neutral), cosmic-ray ionization 
(£=1.3xl0 -17 s -1 ), cosmic-ray induced photo¬ 
dissociation of ice, photolysis, thermal desorption, 
non-thermal desorption (i.e cosmic-ray evapora¬ 
tion), ga s phase and grain surface photodissocia¬ 


tion ( see lHasegawa fe Herbstlll993t iGarrod et al 


120081 : iFogel et ah 2011 : Aikawa et al. 2012 . for 
further details) . I ncident ally, in our analysis, fol¬ 
lowing [Garrodj^llerbst ( 20061) , we assume a CO 
binding energy set to 1150 K (which corresponds 
to a desorption temperature of abo ut ~ 2 1 25K, 
see Supplementary Materials of ICleeves et al. 
20141) . Therefore, the CO snow-line from our 


model is sh i fted by sever al AU in comparison with 
Schwarz fe Berginl ( 20141 ). Although, the chemical 


evolution of the disk is run independently (i.e. 1 D 
vertical model), the code was run for a given range 
of radii resulting in a pseudo two-dimensiona l re¬ 
sult as desc r ibed b y lSchwarz fe Bergin ( 20141) and 
Fog el et al. ( 201ll) . The gas-grain chemical model 
is run for 3xl0 6 years (typical disk age) with 45 
vertical zones and 22 radial zones. 


1 To process faster the data GNU parallel llTangdl2011 ') was 
used along with our code. 


3.1. Disk chemical abundances 


In Figure [TJ we illustrate where the CO snow 
line is in our model by showing the relative 
abundance to molecular hydrogen for CO, COi ce 
and N 2 H+. The latter ion has been reported 
to be a good tracer of the CO snow line since 
it is destroy ed by gas phase CO through reac 


tion m [u l iHersant et al.l 120091 : iQi et al 
Schwarz fc BergiiJ 2014). 


20131 


N 2 H+ + CO ->■ HCO+ + N 2 


( 6 ) 


Our model is consistent with this interpretation. 

In Figure [1] we also show the relative abun¬ 
dance to molecular hydrogen for DCO + along with 
the one of its 2 potential key precursors, CH 2 D + 
and H 2 D+ (see Reactions [2] and [3]). It is immedi¬ 
ately apparent that H 2 D + probes the CO snow 
line as expected from reaction [2] Indeed, the 
H 2 D+ distribution is clearly anticorrelated with 
the one of CO since H 2 D+ is located beyond the 
CO snow line (i.e where CO is adsorbed on ice 
grain mantles) and in the upper disk layers where 
the CO is photodissociated. As for CH 2 D+, it 
is more abundant in the upper disk layers where 
ionization supplied by X-rays is higher. Regard¬ 
ing the overall distribution of the DCO + abun¬ 
dance, it is more or less correlated with that of 
CO, with one important exception: there is a 
greater concentration in the surface of the inner 
disk. More specifically, the DCO + distribution 
is not restricted to the CO snow line but rather 
spreads up to the disk surface with a maximal 
abundance reached in the inner warmer (T~ 50- 
100K) surface layers (see Fig. [1]). 


3.2. On the DCO+ formation pathways 

Reaction [21 involving H 2 D + , is usually consid¬ 
ered as the main s ource of DCO + in dis ks at low 
temperatures (see Mathews et al. 1 120131) . If this 
is the case, the qualitative distribution of DCO + 
abundance should resemble that of CO. However, 
as shown in FigurelU the bulk of the DCO + abun¬ 
dance is located in a narrow layer near the sur¬ 
face. This finding suggests that reaction [2] is not 
the only formation route for DCO + in disks. This 
is further illustrated in Figure [2] that shows the 
DCO + , HCO + , CO, oe and H 2 D + column densities 
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as a function of the radius when reaction [2] is acti¬ 
vated (top panel) and turned off (bottom paneljl. 
The DCO + profile remains almost the same in 
the inner 60 AU whether the H 2 D + channel (reac¬ 
tion [2J is involved or not. Further out, reaction [2] 
contributes for about 40 % in the DCO + abun¬ 
dance, the other main contribution being through 
HCO + (via the reaction HCO + + D —>• DCO + + H, 
Aikawa et al. in prep.). 

One notable feature of Figures 0 and Q] is that 
the formation of DCO + in the inner surface lay¬ 
ers inside of 60 AU involves CH 2 D + or one of its 
product that have previously been created by re¬ 
action [31 Indeed, Figure [3] shows the abundance 
of the DCO + ion for up to the 60 AU when re¬ 
action [3] is activated (left panel) and turned off 
(right panel). It is immediately apparent that re¬ 
action^ which involves the CH 2 D + channel, leads 
to the observed maximal DCO + abundance in the 
inner surface layers of the disk. Likewise, a sim¬ 
ilar result is observed on the abundance profiles 
of DCN and HDCO (see Figure [3]), molecules that 
can both also be formed through reaction[3] In ad¬ 
dition, in Figure U we give the radial distribution 
of some key ions (such as CH 2 D + and CH 4 D+) 
that can lead to the formation of DCO + along 
with the one for H 2 D + when reaction [3] is acti¬ 
vated (top panel) and turned off (bottom panel). 
A decrease in the DCO + , CH 2 D + , dUD -1- and 
CH 4 D + column density is clearly observed when 
reaction[3]is not activated in our network. Finally, 
we find that more than 80% of the DCO + abun¬ 
dance in the inner 60 AU surface layers is due to 
a DCO + formation through the reactions [31 [3] and 
[5] pathway. 

3.3. DCO+ abundance and exothermicity 

In Figure [5l we show that the use of the new 
exothermicities leads the warm channel (i.e. re¬ 
action (3[) to dominate in the DCO + production. 
More specifically, Fig. [5] illustrates the impact, 
within the inner disk, of the change of exother¬ 
micity on the DCO + /HCO + abundance ratio but 
also on that of other easily observable molecules 
(DCN, formaldehyde). One notable feature is that 
the abundance of DCO + , DCN and HDCO is en- 


2 When switching off reaction [2] we keep reactions Hjj" + 
CO -> HCO+ + H 2 and H 2 D+ + CO —► HCO+ + HD 
unchanged. 


hanced in the 20 AU inner radii when the chemical 
modelin g incl udes the new exothermicities set by 
Roueff et al. ( 2013 . AE of 654 K for reaction [3]) as 
opposed to the chemical modeling that incl udes 
the p revious ones (AE = 370 K, see ISmith et al. 
19821) . Furthermore, it is immediately apparent 


that within the 20 AU inner radii the distribution 
of the DCO + /HCO + abundance ratio is affected 
by the change of exothermicity: the latter does 
not increase with increasing r adius w hile usin g th e 


most recent exothermicities ([Roueff et al . 20131) . 


wherea s assuming a lower exothermicity (i.e. the 


one of Smith et al. 19821) the ratio shows a rise 


with increasing radius in the inner disk. 


4. Discussion 

4.1. The effect of the exothermicity in the 
CH 2 D + formation 

Our results are a consequence of the increase 
in the ex otherm icity A E of reaction [31 as sug¬ 
gested by [Roueff et al.1 ( 2013 ). which favors deu¬ 
terium fractionation in warm gas, up to tempera¬ 
tures even as high as T < 300 K. We note that in 
the case of DCO + only, which has its main forma¬ 
tion route from channel [5] and, for an exothermic¬ 
ity of 654 K, a temperature >71 K is hot enough 
to turn off the DCO + enrichment toward the inner 
disk. Indeed, the activation of the backward re¬ 
action is shifted to higher temperatures, enabling 
efficient deuteration closer to the star. This is a 
nice illustration of the influence of the energy bud¬ 
get of deuteration reactions. This is similar to the 
well known effect of the H 2 ortho/para ratio o n 
deuterium chemistry (see e.g. Flower et al. 200(jl : 
at low temperature, the backward counterpart of 
reaction |T| is strongly favored for ortho-H 2 over 
para-H 2 , due to a lower AE for the former. How¬ 
ever, we did not consider the spin states of H 2 in 
the present study. DCO + being mainly produced 
by reaction [31 , we do not expect this omission to 
be of crucial importance in the present case (see 
Sect. 4.4 for further details). 


4.2. 


DCO + a probe of ionization rather 
than the CO snow line 


DCO + has been pr oposed to be a good tr acer 
of the CO snow line bv lMathews et ah . (2013), as¬ 
suming H 2 D+ as its parent molecule. However, 
our study points out that the reaction which in- 
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volves H 2 D+ as a parent molecule (see reaction [2j) 
is unlikely to be the main formation route that 
leads to the bulk of the DCO + abundance within 
the disk. Alternatively, our study shows that re¬ 
actions involving CH 2 D + as a parent molecule of 
DCO + (see reaction [3] and reactions [4] and [5] as 
possible pathways for the formation of DCO + ) 
are responsible for the observed high abundance 
of DCO+ in the inner surface layers of the disk. 
This finding leads us to argue that, instead of tem¬ 
perature, DCO + is rather tracing ionization. In¬ 
deed, the distribution of the DCO + abundance is 
similar to the one of CO except for its high abun¬ 
dance in the the inner warmer surface layers of the 
disk where ionization through X-rays is present. 
In that respect, we predict that DCO + is a tracer 
of the inner tens of AU of protoplanetary disks 
surrounding X-ray active T Tauri stars. 


4.3. Implications: on the production of 
other deuterated molecules 

Chemical pathways involving the CH 2 D + chan¬ 
nel (reaction [3]) can sufficiently lead to the forma¬ 
tion of different deuterated molecules, such as 
DCO + , DCN, HDCO in warm con ditions (e.g . 


T~70-100 K, see Woottenl Il987t lOberg etal 


2012( 1. In that context, our study not only shows 
that the high abundance of DCO + in the inner 
disk results from pathways involving the CH 2 D + 
ion as a parent molecule, but in addition shows 
that the DCN and HDCO abundance patterns 
are also affected by the change of exothermicity 
(see Figs [3] and©. Indeed, the DCN and HDCO 
abundances are enhanced while using the CH 2 D + 
formation route with the most recent exothermic- 


ities (AE=654 K. IRoucff ct al . 20131 ). 


4.4. Implications: on the specific struc¬ 
ture of protoplanetary disks 

Interesti ngly enough, fr o m th eir ALMA ob 


servations, Mathews et al.l ( 2013h analyzed the 
DCO + emission in the disk surrounding HD 163296 
as coming from a narrow (confined between 110 
and 160 AU) ring-like structure. Although the ra¬ 
dial extent of the ring may be subject to revision 
(Qi et al., in prep), such a narrow ring, located 
far from the star, is difficult to reconcile with our 
modeling . Likewise, based on SMA observations 
( Qi et al.ll2008l) of the T W Hydrae di s k and em¬ 
pirical fitting of the data, Oberg et al. ( 2012t ) find 


a central depression and inferred a DCO + abun¬ 
dance profile that increases with abundance with 
the radius. These suggestions disagree with our 
results, si nce we find that wit h a higher activation 
barriers dRoueff et al.l l2013h the warm channel 
(reaction© dominates, leading to a strong abun¬ 
dance of DCO + in the inner regions of the disk. 


On e possible cav e at c ould be the ortho:para 
ratio. Teague et al. ( 2015h have investigated the 
effect of the ortho:para ratio on the deuterium 
fractionation and have found that, by a million 
years, the o:p ratio is low. However, assuming a 
higher o:p ratio might change our r esults by re¬ 
ducing the deuterium enhancement. iRoueff et al.l 
( 20131 ) show that the exothermicity in reaction [3] 
is reduced if there is significant amount of orth o- 


H 2 


present in the gas (see Table 3 of lRoueff et al.l 


20131 ). As a result, the backward reaction is ac¬ 


tually comparable to the old warm channel rates 
(ISmith et al.lll982lf . Thus, one way to chemically 
mimic a hole in the DCO + distribution is to have a 
gradient in the o:p ratio in the disk which is higher 
in the inner regions and lower in the outer disk. 


Another possible caveat could be the ionization 
structure. The presence of structure in the inner 
disk could perhaps prevent the propagation of the 
X-ray photons in the inner disk and may lead to 
a reduced production of DCO + . However, this 
effect may not be as important for the more sub¬ 
stantially flared outer disk. 


Finally, the disk surrounding HD 163296 is 
much warmer than our TW Hydrae disk model. 
In this disk, the CO sn ow line is located around 
150 AU ( Qi et al. 2011 1. while it is around 20 AU 
in our model. Another possibility would thus 
be that the HD 163296 disk is warm enough 
(T>71 K) and/or present a higher o:p ratio to en¬ 
able endothermic reactions to prevent deuterium 
fractionation in the inner 100 AU, thereby trun¬ 
cating the DCO + layer. Direct modeling of the 
DCO + chemistry in the warm disk surrounding 
HD 163296 would be necessary for further com¬ 
parisons. In any case, further observations of 
HD 163296 and TW Hya are required to under¬ 
stand the DCO+ distribution, at the light of our 
new DCO + formation pathway. 
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5. Conclusions 

To summarize, our study shows that the high 
abundance of DCO + in the inner surface layers of 
the disk results from CH 2 D+ as a parent molecule. 
The latter being efficiently formed through reac¬ 
tion [3] in warm conditions (typically T<71 K). As 
a consequence we suggest that DCO + is a reliable 
tracer of ionization in the inner warm disk surface. 
Furthermore, the change in exothermicity, in the 
instance of a low o/p ratio, leads to enhancements 
of other tracers including DCN and HDCO. 

This work was supported by the National 
Science Foundation under grant 1008800. FH 
thanks the Department of Astronomy (University 
of Michigan) for its hospitality. 
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Fig. 1.— Relative abundance to molecular hydrogen. From top left to bottom right: CO, COi ce , DCO" 1 ", 
CH 2 D+, H 2 D+ and N 2 H+. 
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Fig. 2.— Radial distribution of molecular column densities with the H 2 D + channel (reaction [21 top panel) 
and without the H 2 D+ channel (bottom panel). Top panel: The CO; ce , DCO + , HCO + and H 2 D + distribution 
are indicated in dark solid lines, red dashed lines, blue dash dot dot lines and green dotted lines, respectively. 
Note that in both panels the COi ce column density is divided by a factor 10 5 . 
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Fig. 3. DCO + (top), DCN (middle) and HDCO (bottom) abundance relative to molecular 
panels: with the CH 2 D + channel (reaction [3]) . Right panels: without the CH 2 D + channel. 
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Fig. 4.— Radial distribution of molecular column densities with the CH 2 D + channel (reaction [3] top panel) 
and without the CH2D + channel (bottom panel). Top panel: The COi ce , DCO + and HCO + distribution are 
indicated in dark solid lines, red dashed lines and blue dash dot dot lines, respectively Bottom panel: The 
CH 3 D+, H 2 D + , CH4D+ and CH 2 D + are indicated in yellow solid lines, green dashed lines, magenta dotted 
lines and cyan dash dot dot lines, respectively. Note that in both panels the COj ce column density is divided 
by a factor 10 5 . 
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Fig. 5.— Impact on the exothermicity on the DCO + /HCO + (dark line), DCN/HCN (cyan dashed line) 
and HDCO/H 2 CO (blue dotted line) column density ratios. Top panel: Abundance ratio measur e d from 
modeling that includes the new exothermocity of reaction [3] (AE = 654 K) set bv iRoueff et al.l (2013:). 
Bottom pa nel: Abunda nce r atio measured from modeling that includes the previous exothermicity (AE = 
370 K, see lSmith et al. 1982). 
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